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Abstract. The composition of ultra-liigii energy 
(UHE) cosmic rays E > 10^"^ eV is still unknown. 
The observation of UHE photons would extend the 
observed electromagnetic spectrum to highest energy 
and open a new channel for multimessenger obser- 
vations in the universe. Current limits on the photon 
flux already constrain "exotic" scenarios where a 
large number of photons is expected by the decay 
products of supermassive X-particles. Motivated by 
the growing exposure of UHE cosmic ray experiments 
- like the Pierre Auger Observatory - the observation 
of conventionally produced GZK photons may be in 
reach in the near future. We investigate UHE particle 
propagation using the Monte Carlo code CRPropa. 
Particularly, the expected photon fluxes normalized 
to current experiments as weU as prospects for 
future experiments are Ulustrated. Varying source 
and propagation scenarios are analyzed and the 
impact on secondary GZK photons is shown. For 
the specific case of Centaurus A, we study which 
source parameters can be tested by searching for the 
expected GZK photons. 

Keywords: cosmic ray propagation, UHE photon 
flux, Centaurus A 

I. Introduction 

The origin and nature of the highest energy cosmic 
rays (E > 10^^ eV) is still one of the most pressing 
questions of astroparticle physics. However recent de- 
velopments show a clear evidence of a suppression in 
the cosmic ray flux at highest energies. HiRes reported 
the observation of the GZK cutoff above ^ 6 • 10^^ eV 
with 5 standard deviation significance fT |. Furthermore, 
the Pierre Auger Observatory rejects the hypothesis that 
the cosmic ray spectrum continues with a constant slope 
above 4 • 10^^ eV, with a significance of 6 standard 
deviations [2J. 

The composition at these energies still remains a 
mystery. The Pierre Auger Observatory revealed a corre- 
lation between the arrival directions of ultra high energy 
cosmic rays (UHECR) with energy above 6 • 10^^ eV 
and the positions of active galactic nuclei (AGN) within 
^ 75 Mpc L3J. This perhaps indicates a lighter composi- 
tion since heavier nuclei are more effected by magnetic 
fields. However, measurements of the depth of shower 
maximum Xj^ax of air showers seem to indicate also a 
heavier component H. 

In either case, energy loss by propagation effects limit 




energy [eV] 

Fig. 1. Spectrum of secondary pliotons generated by pion and pair 
production from a single UHECR proton source at a given distance. 
We consider here a one-dimensional model, with an injection spectral 
index a = 2.5 and maximum energy of lO^^'^ eV. No magnetic fields 
were taken into account. At a source distance of ~ 10 Mpc most of 
the UHE photons are produced. For closer distances the EM cascade's 
development has insufficient time to produce a sufficient number of 
UHE photons whereas for large distances the UHE photon population 
may cascade down to lower energies (see also [,5]). 

the UHECR horizoii] distance to below ^ 70 Mpc at 
energies > 10^° eV and give rise to secondary particle 
production. 

To get a clue of an answer of the raised questions 
it is therefore desirable to expand the knowledge of 
particle propagation through the local universe. The 
photon background is a key ingredient for understand- 
ing the properties of particle propagation. At energies 
> 5-10^^ eV the main channel of energy loss for primary 
protons is photo-pion production in interactions with 
background radiation fields which generates the already 
mentioned GZK feature. Here, the low energy photon 
can Lorentz transform into a 7-ray in the rest frame of 
a very-high energy particle. The cross section increases 
strongly at the A+(1232) resonance. The process can 
be described as 

p + 7^A+(1232) ^ n + TT+ 

p + tt" . 

In addition, also further baryon resonances can be 
excited at increasing energy. The produced neutral pions 
decay into two UHE photons which in turn are distance 
limited by 77 interactions with background photons. 

In Fig. [T]a proton source with spectral index a — 2.5 
was simulated at various distances. The simulations 
were made using the numerical tool CRPropa [61 which 

'Here the horizon d is defined as the distance within which 90% 
of arriving particles originated. 
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Fig. 2. Simulation matrix of varying source parameters. A changing spectral index is shown on the y-axis in combination with a varying 
maximum energy of the source. Color coded is the expected photon flux F^OEcV above a threshold energy of 10 EeV. The upper limit on the 
flux of photons above 10 EeV derived in j8| is i^i'nS'V ~ "^'^ ' km~^ sr~^ yr~^. The white dashed line indicates a photon flux level 



of i^i''rIS!,Y/30, where the factor 30 is a rough estimate of how much the final total exposure of Auger South exceeds the exposure used for 



XOEeV' 
mit 
lOEeV 



^InS'A^ ■ One sees that interesting combinations of a and i?max can be tested by searching for the expected photon flux. 



is described in more detail in Sec. |ll] The resulting 
UHE photon flux by pion and pair production processes 
is shown. A region around the source exists where 
the UHE photon flux is maximal. For closer distances 
(e.g. 2 Mpc) the GZK effect does not yet efficiently 
produce UHE photons, whereas for larger distances 
(e.g. 50 Mpc) the UHE photon population may go 
into the development of a full electromagnetic cascade 
with the main flux arriving at GeV-TeV energies. UHE 
photons can therefore provide information on local 
UHECR sources. 

Up to now no UHE photons were observed yet |7|. 
The strongest constrains on the UHE photon flux at 
£^ > 10 EeV were set by the Pierre Auger Observatory 
using surface detector data lIS). Recently, new photon 
limits down to 2 EeV were presented using a combi- 
nation of surface and fluorescence detector |9|. These 
results are statistically limited and can be improved in 
the future with the capability of measuring the diffuse 
flux as well as identifying discrete sources. It has also 
been suggested that a major UHECR flux may arise from 
just a few nearby AGN such as Centaurus A ITOll . At 
a distance of 3.4 Mpc HI J Centaurus A is by far the 
nearest active radio Galaxy. 

II. CRPropa simulation 

The interplay between different astroparticle physics 
experiments has become very important. Existing and 
planned projects range from UHECR observations 
like the Pierre Auger Observatory, to neutrino 
telescopes iTTSI . lfT3l . as well as ground and space based 
7-ray detectors operating at TeV and GeV energies, 
respectively lfT4l . Even if a putative source were to 
produce exclusively UHECR, photo-pion and pair 
production by protons on the photon background would 
lead to guaranteed secondary photon and neutrino 
fluxes that could be detectable. With this motivation a 



numerical tool called CRPropa Q has been developed 
that can treat the interface between UHECR, 7-ray and 
neutrino astrophysics, and large scale magnetic fields. 

Pion production is modelled by using the event gen- 
erator SOPHIA 1 15 1 that has been explicitly designed to 
study this phenomenon and is augmented in CRPropa for 
interactions with a low energy extra-galactic background 
light (EBL). Unlike pion production, pair production by 
protons is taken into account as a continuous energy loss 
due to the low inelasticity. 

The EM cascade code is based on ||16i . All relevant 
interactions with background photons are taken into ac- 
count and implemented in CRPropa including single pair 
production, double pair production, inverse Compton 
scattering and triplet pair production. There are three 
different photon backgrounds implemented in CRPropa. 
The most important one is the cosmic microwave back- 
ground (CMB) with a well known redshift evolution. 
For the infrared background a model of Primack et al. 
is used |17|. This becomes important for EM cascades 
around the threshold for pair production and is less 
significant in the UHE region. Above ~ 10^^ eV interac- 
tions with the universal radio background (URB) become 
more important where it can inhibit cascade development 
due to the resulting small pair production length. We use 
a model based on observations IJSil. 



In this study we investigate one-dimensional particle 
propagation for distances of 3.4 Mpc i.e. the distance 
to Centaurus A. Here we assume a proton source with 
spectral index a and maximum energy i^max- Protons 
are injected and propagated (assuming no magnetic 
field) towards the observer The resulting EM spectra 
are recorded and weighted according to the Auger flux 
spectrum 1 19| as follows: we assume in a first simplified 
step that the total flux above 57 EeV originates from 
a source at a distance of 3.4 Mpc, i.e. in this case 
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Fig. 3. Simulation matrix of varying source parameters. A changing spectral index is shown on the j/-axis in combination with a varying 
maximum energy of the source on the x-axis. Color coded are the expected number of secondary photons above a threshold energy of 10 EeV 
for an exposure of 3130 km^ sr yr, corresponding to | 8 |. The white dashed line indicates one photon above 10 EeV. The red numbers indicate 
an increased exposure by a factor of 30, coiresponding to the expected sensitivity of Auger South. 



27 hadron events above 57 EeV with 7000 km^ sr yr 
exposure (as observed in [,19|). The ratio /?7ooo of the 
observed 27 events to the simulated total number of 
arriving particles above 57 EeV is then also used to scale 
the simulated photon flux above a certain energy E (here 
we take E = 10 EeV), i.e. 



ArlOEcV _ atIOEcV 



7000 



(1) 



where N}^'^^^ is the simulated number of photons above 
10 EeV and A^.^'^cra '^he number of photons that are 
expected to be observed above 10 EeV with an exposure 
of 7000 km2 sr yr [19 1. 

The expected integrated 7-flux above 10 EeV, F]^'^^'^^ , 
is then calculated via 



uilOEcV 
-^7 



ArlOEeV 
^'^7,7000 



7000 km^ sr yr 
III. Results 



(2) 



In Fig. I2] the expected photon flux p^^EcV gjjown 
for varying source parameters a and i^max- Compared 
to the upper limit on the flux of photons above 10 EeV 
derived in 18J of i^|'J5^'*v = 3.8-10"^ km^^ si--i yi--i ah 
simulated source parameter combinations are compatible 
with the current upper limit. That is, the present upper 
limit on the photon flux does not yet constrain Centaurus 
A as a strong source of UHE protons. We checked 
that a constant tans verse magnetic field of 100 pG has 
just a marginal effect on the UHE photon flux in this 
scenario. The white dashed line in Fig. |2] indicates an 
improved photon flux level of Fl'^^^^ /30, corresponding 
to a rough estimate of the expected sensitivity of Auger 
South. Here, certain parameter combinations of a and 
^-max produce a larger photon flux and can thus be 
tested. 

The number of arriving photons above 10 EeV for 
varying source parameters that are expected to be ob- 
served with an exposure of 3130 km^ sr yr (correspond- 



ing to the exposure used for F^^^y H) 

^7!fi30 = ^7"'"=^ • 3130 km^ sr yr (3) 

is shown in Fig. [3] For large i^max and small a the 
expected number increases up to a few photons. For 
an increased exposure of a factor of 30 the number of 
arriving photons also scales by a factor of 30 as indicated 
by the red numbers in Fig. [5] One sees that, depending 
on the source parameters, up to several 10 events could 
be expected. 

So far, we made the simplified assumption that the 
total flux above 57 EeV is produced by the source — 
irregardless of the question whether the observed shape 
of the spectrum or the distribution of arrival directions 
were reproduced. Relaxing now this assumption one can 
ask for the fraction 



l^limit / V 
^lOEcV/^ 

FWEeV^a,E^,^) 



(4) 



where X — 1 refers to the exposure used for -fjoEcV 
and X ~ 30 to the Auger South sensitivity. Assuming 
now an observed photon upper limit Fl]^^^^y/X, for / > 
1, no constrain on source parameters is possible, while 
values of / < 1 indicate which fraction of the total 
flux would still be allowed from the source for a given 
combination of a and E'max- As is clear from Fig. |2l 
/ > 1 for X = 1, for all simulated combinations of a 
and £;,nax- 

In Fig. m the case of X — 30 is shown, and the 
corresponding fractions of the total flux still allowed by 
the source can be extracted. For instance, in case of 
a ~ 2 and i?,„ax — 10^^ eV, no more than ^ 30% of 
the total cosmic ray flux could be due to protons from 
Centaurus A. 

IV. Summary 

We studied GZK photon fluxes expected at Earth 
using the MC program CRPropa. Due to the competition 
between production of GZK photons (increasing with 
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Fig. 4. Simulation matrix of varying source parameters. Color coded is the fraction / : 
of Fig. [2). 



(^]oEcv/30)/i^^''"^"^(a, En^ax) (see also caption 



travel distance of the mother nucleon) and attenuation 
of GZK photons (increasing with travel distance), the 
expected photon fluxes show a non-trivial dependence 
on the source distance (Fig. [T] see also e.g. Il20ll ). 

Regarding the specific case of Centaurus A (see 
also ||2TI ). the current photon flux limit ^ does not 
yet constrain Centaurus A as a strong source of UHE 
protons for the investigated range of spectral indices 
a and maximum energies i?max (Fig. EJ- However, the 
sensitivity that will be accumulated by Auger South will 
allow interesting constrains for a broad range of a and 
i^max (Fig. El and Fig. HI. 

Depending on source parameters, the number of GZK 
photons above 10 EeV may reach several 10 over the 
lifetime of Auger South (Fig. O. We conclude that the 
search for UHE photons helps to provide significant 
clues about the characteristics of potential astrophysical 
sources. 

We note that we regarded only GZK photons, i.e. 
photons produced during the propagation. The photon 
fluxes may be enhanced in case of interactions at the 
source. Such studies, as well as 3d-simulations with 
CRPropa, are in progress. 
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